arXiv:l504.05814v2 [math.DG] 13 Jan 2016 


A combinatorial Yamabe problem on two and 
three dimensional manifolds 

Huabin Ge, Xu Xu 
January 14, 2016 


Abstract 

In this paper, we define a new discrete curvature on two and three dimensional 
triangulated manifolds, which is a modification of the well-known discrete curvature 
on these manifolds. The new definition is more natural and respects the scaling 
exactly the same way as Gauss curvature does. Moreover, the new discrete curvature 
can be used to approximate the Gauss curvature on surfaces. Then we study the 
corresponding constant curvature problem, which is called the combinatorial Yamabe 
problem, by the corresponding combinatorial versions of Ricci flow and Calabi flow 
for surfaces and Yamabe flow for 3-dimensional manifolds. The basic tools are the 
discrete maximal principle and variational principle. 
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1 Introduction 

1.1 Background and Preliminaries 

One of the central topics in differential geometry is the existence of canonical metrics on a 
given smooth manifold, especially metrics with constant curvatures. One can ask similar 
questions on triangulated manifolds. Different from the standard processes on smooth 
manifolds, where all calculations are done in local coordinates, Regge [H] suggested an 
entirely different way, where the calculations are done in a geometric simplex. Since 
the geometric information in a single simplex is included in the lengths of edges, the 
geometric approach to a triangulated manifold may be started with a piecewise flat metric, 
which assigns each edge a Euclidean distance such that each simplex can be realized as 
a Euclidean simplex. Piecewise flat metric brings singularities on codimension two sub- 
simplices, which are explained as simplicial curvature tensors in p El Eg [39]. Besides 
defining edge lengths directly, one can also define discrete metrics at vertices and derive 
the edge lengths indirectly. In fact, Thurston [ig introduced the notion of circle packing 
metric on a triangulated surface, while Cooper and Rivin introduced the notion of sphere 
packing metric on a triangulated 3-manifold. We shall study the circle (sphere) packing 
metrics and the corresponding discrete curvatures in this paper, and we want to know 
when there exist canonical discrete metrics and how to find them. 

Suppose M is a closed surface with a triangulation T = {V,E,F}, where V,E,F 
represent the sets of vertices, edges and faces respectively. Let <I> : E —)> [0, |] be a function 
assigning each edge {ij} a weight ^ij G [0, |]. The triple (M,T, 4*) will be referred to as 
a weighted triangulation of M in the following. All the vertices are ordered one by one, 
marked by ui, • • • , vm, where N = V'^ is the number of vertices. We use z ~ j to denote 
that the vertices i and j are adjacent if there is an edge {ij} G E with i, j as end points. 
Throughout this paper, all functions / : R —)• IR will be regarded as column vectors in 
and /* is the value of / at i. And we use C{V) to denote the set of functions defined 
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on V. Each map r : E —>■ (0, +oo) is called a circle packing metric. Given (M, T, <h), we 
attach each edge {ij} a length 

lij = ^rf + r‘^ + ^nrj cos ^ij. (1.1) 

Thurston proved [H] that the lengths {lij,ljk,kk} satisfy the triangle inequalities, which 
ensures that the face {ijk} could be realized as an Euclidean triangle AviVjVk with lengths 
{lij,ljk,hk}- In this sense, the triangulated surface could be taken as gluing 

many Euclidean triangles coherently. However, this gluing procedure produces singulari¬ 
ties at the vertices, which are described by the discrete curvatures. Suppose is the inner 
angle of the triangle AviVjVk at the vertex i, the classical well-known discrete Gaussian 
curvature at i is defined as 

K, = 27r- Y. (1-2) 

AviVjVj^GF 

where the sum is taken over all the triangles with Vi as one of its vertices. Discrete Gaussian 
curvature Ki satisfies the following discrete version of Gauss-Bonnet formula |14j : 

YKi = 27rx{M). (1.3) 

i&V 

Denote Kav = it is natural to ask the following combinatorial Yamabe problem: 

Question 1. Given a weighted triangulated surface (M, 7~, 4>), does there exist a circle 
packing metric r that determines a constant i^-curvature? How to find it if it exists? 

Different from the smooth case, the constant iiT-curvature metric, i.e. a metric r 
with Ki = Kav for all z G V, does not always exist. Thurstion [l8] found that, besides 
the topological structure, the combinatorial structure plays an essential role for the 
existence of constant E'-curvature metric. In fact, for any proper subset I C V, let Fj be 
the subcomplex whose vertices are in I and let Lk{I) be the set of pairs (e, v) of an edge 
e and a vertex v satisfying the following three conditioins: (1) The end points of e are not 
in /; (2) u is in /; (3) e and v form a triangle. Thurston [38] proved that 

Theorem 1.1. (Thurston) Given a weighted triangulated surface (M, 7~, <I>), the ex¬ 
istence of constant X-curvature metric is equivalent to the following combinatorial- 
topological conditions 

2vrx(M)^>- (7r-ch(e)) + 27rx(E,), VI : </. g / C y. ( 1 . 4 ) 

' ' {e,v)eLk(l) 

Moreover, the constant H-curvature metric is unique, if exists, up to the scaling of r. 
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In fact, the following result is proved for the classical discrete curvature K IMIEZIIII]. 

Theorem 1.2. (Thurston-Marden&Rodin-Chow&Luo) Given a weighted triangu¬ 
lated surface (M, T, <h), denote 




where 


^GB = |x G 


R 


N 

2 = 1 


Xi = 2 tt 


(1.5) 


( 1 . 6 ) 


and for any nonempty proper subset I C V, 

^/ = |x G IR^I ^ {tt - ^{e))+ 27rx{Fi)y (1.7) 

2G/ {e,v)^Lk{I) 

The space of all admissible iiT-curvatures \^K = K{r)\r G lIl>o} is exactly '3^. 

Chow and Luo |14) first established an intrinsic connection between Thurston’s circle 
packing metric and the surface Ricci flow. They introduced a combinatorial Ricci flow 

-r- = —F-iVi 

dt 

with the normalization 

^ = {Ka, - Ki)n. (1.8) 

Then they reproved Theorem |1.1| by the curvature flow methods. In fact, they proved: 

Theorem 1.3. (Chow&Luo) Given a weighted triangulated surface (M, T, $), the so¬ 
lution to the normalized combinatorial Ricci flow (1.8) converges (exponentially fast to a 
constant iL-curvature metric) if and only if there exists a constant iL-curvature metric. 

Inspired by [TTj, the first author introduced a combinatorial Calabi flow 

Tt = 

in |2nj and proved similar results: 

Theorem 1.4. (Ge) Given a weighted triangulated surface (M, T, <h), the solution to the 
combinatorial Calabi flow (1.9) converges (exponentially fast to a constant IL-curvature 
metric) if and only if there exists a constant iL-curvature metric. 

Theorem |l.lt Theorem 1.2 Theorem 1.3 and Theorem 1.4 answer Question [^perfectly. 
The first two theorems give a necessary and sufficient combinatorial-topological condition 
for the existence of a constant RT-curvature metric, while the last two theorems give two 
different efficient ways to hnd a constant RT-curvature metric, i.e. for any initial value 
r(0), the solutions to the combinatorial Ricci flow and Calabi flow converge automatically 
to a constant RT-curvature metric exponentially fast. 
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1.2 Main results 


It seems that the classical discrete Gauss curvature K is a suitable analogue for smooth 
Gauss curvature. However, there are some intrinsic disadvantages for the classical discrete 
Gauss curvature. For one thing, the smooth Gauss curvature changes under the scaling of 
Riemannian metrics while the discrete Gauss curvature dose not. For another, the discrete 
Gauss curvature K can not be used directly to approximate the smooth Gauss curvature 
when the triangulation is finer and finer. Motivated by the two disadvantages, we modify 
the classical discrete Gauss curvature Ki to 

R^ = ^ ( 1 . 10 ) 

G 

at the vertex i, which is called as i?-curvature in the following. In some sense, this new 
definition is more natural and respects the scaling exactly the same way as the Gauss cur¬ 
vature does. Furthermore, this new definition can be used directly to approximate smooth 
Gauss curvature (see Section [^. It is natural to consider the following combinatorial 
Yamabe problem for R-curvature. 

Question 2. Given a weighted triangulated surface does there exist a circle 

packing metric r that determines a constant i?-curvature? How to find it if it exists? 


We introduce the corresponding combinatorial versions of the Ricci flow and the Calabi 
flow for surfaces and the Yamabe flow for 3-dimensional manifolds using the new definition 
of discrete curvature. We give an answer to Questionby discrete curvature flow methods 
and variational methods. We just state the main results in this subsection. 

We introduce a new combinatorial Ricci flow 

^ = {Rav - Ri)9i, (1.11) 

where Rav = 27rx(M)/||r||2 and gi = r? is a discrete version of Riemann metric ten¬ 
sor. This flow has the same form as combinatorial Yamabe flow in two dimension. By 
introducing a new discrete Laplace operator 


A/. 




where / is a function defined on all vertices and uj = In gj is a coordinate transformation, 
we find that the flow (1.11) exhibits similar properties to the smooth Ricci flow on surfaces. 
For example, the curvature Ri evolves according to a heat-type equation 

= Ai?j -|- Ri{Ri — Rav), 
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which has almost the same form as the evolution of the scalar curvature along the classical 
Ricci flow and permits a discrete maximum principle. 

We can also introduce a new combinatorial Calabi flow 


dgi 

dt 




( 1 . 12 ) 


which is similar to the smooth Calabi flow on surfaces. By considering a discrete quadratic 
energy functional C{r) = where ipi = Ki — RavV^, we can define a modihed Calabi 

flow 

(1.13) 


u = -^VuC. 


Combining Corollary |3.19 Theorem 4.2, Theorem 4.6 Theorem 3.24 and Theorem 


3.26, we have the following main theorem. 


Theorem 1.5. Given a weighted triangulated surface (M, T, 4>) with x(M) < 0. Then 
the constant i?-curvature metric is unique (if exists) up to scaling, and the following 
statements are mutually equivalent: 


(1) There exists a constant i?-curvature circle packing metric; 


(2) The solution to the combinatorial Ricci flow (1.11) converges; 


(3) The solution to the combinatorial Calabi flow (1.12) converges; 


(4) The solution to the modified Calabi flow (1.13) converges; 


(5) There exists a circle packing metric r* such that, for any nonempty proper subset I 
of V, 

^*2 

^ (vr - 4>(e))+27rx(F/); (1.14) 

{e,v)£Lk(I) 




(6) The combinatorial-topological conditions ^ H R^q ^ (j) when x{^) < 0; ^ n 

{0} 7 ^ (f) when x{^) = 0 valid, where ‘3^ is defined in (1.5). 

Remark 1. When x{^) > 0; the conclusions in Theorem 1.5 are not true. Example]^ 
shows that the combinatorial Ricci flow for the sphere with the tetrahedron triangulation 
may not converge, while Exampleshows that there are more than one constant curvature 
metric on the sphere with the tetrahedron triangulation. 

Theorem 1.5 and Remark provide a suitable answer to the Question 
As pointed out in Remark it’s surprising that the combinatorial flows are quite dif¬ 
ferent from the smooth flows. The combinatorial structure of the triangulation bring some 
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extra trouble that needs special considerations. However, for the surfaces with x{^) > 0; 
it’s very interesting that we could evolve the flow ( |1.11[ ) very well by substituting the “area 
element” rf by or, more generally, rf with a < 0. In fact, for any a G R, we can 
dehne the a-curvature as 


RqA — 


K, 


(1.15) 


at each vertex i. The idea behind this definition is to consider r“ as a metric (of a 
order). From the viewpoint of Riemannian geometry, a piecewise flat metric is a singular 
Riemannian metric on or M^, which produces conical singularities at all vertices. For 
any a G M, a metric g with conical singularity at a point can be expressed as g{z) = 
locally. Choosing f{z) = — Ina^, then g{z) = \dz°^\‘^. Comparing r" 
with I dz" I, the a-metric may be taken as a discrete analogue of conical metric to some 
extent. By introducing the corresponding a-flows, we get the following result for constant 


a-curvature problem, which is a combination of Theorem 5.4, Theorem 5.6 and Corollary 

ESI 


Theorem 1.6. Given a weighted triangulated surface (M, T, <!*) with ax(M) < 0. Then 
the constant a-curvature metric is unique (if exists) up to scaling, and the following 
statements are mutually equivalent: 

(1) There exists a constant a-curvature circle packing metric; 


(2) The solution to the a-Ricci flow (5.4) converges; 


(3) The solution to the a-Calabi flow (5.7) converges; 


(4) The solution to the modified a-Calabi flow (5.8) converges; 


(5) There exists a circle packing metric r* such that, for any nonempty proper subset I 
of V, 

27^x(M)^^0^ > - (tt - <P(e)) + 2vrx(Fj); (1.16) 

{e,v)&Lk{I) 

(6) The following combinatorial-topological conditions 

(i) n R^q 7 ^ (/), when a > 0 and x(M) < 0; 

(ii) '3^ n R>q 7 ^ (/>, when a < 0 and x(M) > 0; 

(hi) (1, • • • , 1)^ G when ax(M) = 0. 


are valid, where 3^ is dehned in (1.5) 
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Remark 2. Theorem 1.6 contains Theorem 1 1.1 1.3 1.4 1.5 as special cases. Take a = 2, 


then Theorem 1.6 is reduced to Theorem |1.5t Take a = 0, then the equivalence between 
(1) and (5) is in fact Theorem [13 the equivalence between (1), (2) and (4) is in fact 
Theorem o and Theorem 11.41 


Theorem 1.6 can not be generalized to the case of ax{M) > 0 due to Remark]^ so it 
is optimal to some extent. We expect there are triangulations in this case that admits at 
least one constant a-curvature metric, or more ideally, admits a solution to discrete Ricci 
flow that converges to a constant a-curvature metric. By the discrete maximum principle, 
we have the following existence result for non-negative constant a-curvature metric. 


Theorem 1.7. Suppose (M, T, 4>) is a weighted triangulated surface. If there is a metric 
r satisfying Ra,i > 0 for all i G V, and 

- (tt - 4>(e)) + 27rx(T,) < 0, VI : 0 C / C y, 

{e^v)£Lk(I) 


then there exists a non-negative constant a-curvature metric r*. 

For a compact manifold with a triangulation T = {V, E, F,T}, where T is the set of 
all tetrahedrons. A ball packing metric r : R —)• (0, -|-oo) evaluates each {i, j} G E a length 
kj = ri + rj. Denote aijki as the solid angle at a vertex i in the tetrahedron {i,j, k, 1} G T. 
Cooper and Rivin [18] defined a discrete scalar curvature (we call it “CR-curvature” for 
simple in the following ) at each vertex i as 

dvr ^ ^ Oiijki- 


To study the constant CR-curvature problem, Glickenstein |24[ 125] first introduced a 
combinatorial Yamabe flow 


dvi 

dt 


-Kiri. 


(1.17) 


Inspired by Glickenstein’s work, the first author of this paper and Jiang |21) modify Glick- 

uV • 

-^ = {X-Ki)r,, (1.18) 


enstein’s flow (1.17) to 


where A = get better convergence results. Similar to the two dimen¬ 

sional case, we give a new definition of the combinatorial scalar curvature 


Ri = 


K^ 



Similar to the smooth case, we find that the sphere packing metrics with constant com¬ 
binatorial scalar curvature are isolated and are exactly the critical points of a normalized 
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Einstein-Hilbert-Regge functional ■ Hence we propose to study a com¬ 

binatorial Yamabe problem with respect to the new combinatorial scalar curvature Ri. 


Question 3. Given a 3-dimensional manifold M with triangulation T, find a sphere pack¬ 
ing metric with constant combinatorial scalar curvature in the combinatorial conformal 
class SOTt-. Here is the space of admissible sphere packing metrics determined by T. 


Set Rav = X] we introduce the following combinatorial Yamabe flow 


^ = {Rav - Ri)gi 


(1.19) 


to study the combinatorial Yamabe problem. Nonpositive constant curvature metrics are 


local attractors of the flow (1.19), which implies the following result. 


Theorem 1.8. Suppose r* is a sphere packing metric on {M,'T) with nonpositive con¬ 
stant combinatorial scalar curvature. If ||r’(0) — is small enough, the solution of the 
normalized combinatorial Yamabe flow (6.14[) exists for all time and converges to r*. 


The paper is organized as follows. In Section We introduce the new definition of 
discrete Gauss curvature for triangulated surfaces with circle packing metrics. In Section]^ 
We introduce the combinatorial Ricci flow and use it to give a solution to the 2-dimensional 
Yamabe problem and the prescribing curvature problem. In Section we introduce the 
combinatorial Calabi flow and study its properties. In Section we introduce the notion 
of a-curvature and the corresponding a-flows and then study the corresponding constant 
curvature problem and prescribing curvature problem using the a-flows. In Section we 
introduce the new definition of combinatorial scalar curvature on 3-dimensional manifolds 
and give a proof of Theorem 1.8 In Section We list some unsolved problems closely 
related to the paper. 


2 The definition of combinatorial Gauss curvature 

The well-known discrete Gauss curvature Ki has been widely studied in discrete geometry, 
we refer to fiHETlES] for recent progress. However, there are two disadvantages of the 
classical definition of Ki. For one thing, classical discrete Gauss curvature does not perform 
so perfectly in that it is scaling invariant, i.e. if fj = Ar^ for some positive constant A, then 
Ki = Ki, which is different from the transformation of scalar curvature R\g = \~^Rg in 
the smooth case. For another, classical discrete Gauss curvature can not be used directly 
to approximate smooth Gauss curvature when the triangulation is finer and finer. As the 
triangulation of a fixed surface is finer and finer, we can get a sequence of polyhedrons 
which approximate the surface. However, the classical discrete Gauss curvature at each 
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Figure 1: a local triangulation of sphere 


vertex i tends to zero, for all triangles surrounding the vertex i will finally run into the 
tangent plane and the triangulation of the fixed surface locally becomes a triangulation of 
the tangent plane at vertex i. To approximate the smooth Gauss curvature, dividing Ki 
by an “area elment” Ai is necessary. Since we are considering the circle packing metrics, 
we may choose the area of the disk packed at i, i.e. Ai = yrr?, as the “area elment” 
attached to the vertex i. Omitting the coefficient vr, we introduce the following definition 
of discrete curvature on triangulated surfaces. 


Definition 2.1. Given a weighted triangulated surface (M, T, ‘h) with circle packing met¬ 
ric r : F —)• (0, Too), the discrete Gauss curvature at the vertex i is defined to be 


Ki 

Ri = ^, 

n 


where Ki is the classical discrete Gauss curvature defined as the angle deficit at i by (1.2). 


The curvature Ri still locally measures the difference of the triangulated surface from 
the Euclidean plane at the vertex i. In the rest of this subsection, we will give some 
evidences that the Definition |2.1| is a good candidate to avoid the two disadvantages we 
mentioned in the above paragraph. Firstly, let us have a look at the following example. 


Example 1. Gonsider the standard sphere imbedded in The smooth Gauss curva¬ 
ture of is -|-1 everywhere. Gonsider a local triangulation of at the north pole, which 
is denoted as a vertex i. O is the origin of R^. /i is a point lying in Oi, see Figure The 
intersection between the horizonal plane passing through h and is a circle. Divide this 
circle into n equal portions. Denote any two conjoint points by j and fc, thus i, j, k forms 
a triangle. Denote x = AiOj. Taking rj=rj=rfc=ij'/2=sin |, ^ij = = 0 and suitably 

choosing ^jk, we can always get a circle packing metric locally with Ijk = 2 sin ^ sin x. 
Moreover, we can get an infinitesimal triangulation at vertex i by letting n —?• -foo and 
X —)• 0. By a trivial calculation we can get = Zkij = 2 arcsin(sm^cos|), hence the 
classical discrete Gauss curvature is Ki = 2'K — 2narcsin(sm^cos|). It’s obviously that 


lim Ki = 0. 
n^oo, 
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However, the discrete Gauss curvature Ri approaches the smooth Gauss curvature in the 
following way 


Ri _ Ki 

vr vrr? 


27r — 2 n arcsin(sin-cos I) 


TT sin^ I 


+ 1 . 


Secondly, let’s analyse from the view point of the scaling law of Riemann tensor and 
Gauss curvature. Recall that, for a curve 7 : [a, 6 ] —>■ M in a Riemannian manifold 
the length of the curve is defined to be L{'y,g) = g{'y{t),j{t))dt, where j{t) 

is the tangent vector of 'y{t) in M. So we have L{'y,g) = X^^'^L{j,g) g = Xg for some 
positive constant A. Note that, for a weighted triangulated surface (M, T, <h) with circle 
packing metric r : R —)■ (0, + 00 ), the length lij of the edge {ij} is given by (1.1). So the 
quantity corresponding to the Riemannian metric g on the weighted triangulated surface 
(M, T, with circle packing metric r should be quadratic in r, among which r? is the 
simplest. Furthermore, according to the definition Rj = ^, we have Ri{Xrf) = X~^Ri{rf), 
which has the same form as the transformation for the smooth scalar curvature. 


Finally, the Definition 2.1 is geometrically reasonable. Let us recall the original def¬ 
inition of Gauss curvature. We just give a sketch of the definition here and the readers 
could refer to m for details. Suppose M is a surface embedded in and is the unit 
normal vector of M in IR^. v defines the well-known Gauss map v : M ^ <Z R^. Then 
the Gauss curvature at p G M is defined as 

Aera v{A) 


K{p) = lim 

A-^p 


Aera A 


where the limit is taken as the region A around p becomes smaller and smaller. For 
the weighted triangulated surface (M, T, $) with circle packing metric r, assume it is 
embedded in and take the normal vector as a set-valued function, then the definition 
of discrete Gauss curvature R* = ^ is an approximation of the original Gauss curvature 

i 

up to a uniform constant vr, as the numerator Ki is the measure of the set-valued function 
V and the denominator rf is the area of the disk with radius rj up to a uniform constant 


Note that the classical discrete Gauss curvature satisfies the Gauss-Bonnet identity 


(1.3). If we define a discrete measure p on the vertices by pi = rf, then we have 

N 


Jm 


( 2 . 1 ) 


for any / G C{V). Using this discrete measure, we have the following discrete version of 
Gauss-Bonnet formula 

I Rdfi = Rjdpi = 27rx(M). (2.2) 


IM 
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In this sense, the average curvature is 


D _ /m _ 2vrx(M) 

■^av — n 7 — I I I 19 ’ 

Irlr 

where Hrlp = r? is the total measure of M with respect to /r. 

We give the following table for the corresponding quantities for smooth surfaces and 
triangulated surfaces. 



Smooth surfaces 

Weighted triangulated surfaces 

Metric 

g = gijdx^dx^ 

9i = rf 

Length 

Hi, 9) = fa V9{i{t),i{t))dt 

kj = +'^j + cos ^ij 

Measure 

d/i = Y^det 9ijdx 

dm = rf 

Curvature 

Gauss curvature K 


Total curvature 

= 27rx(M) 

J2i Ridm = 2 ttx{m) 


Table 1: Corresponding between smooth surfaces and triangulated surfaces 


For the curvature Ri defined by (2.1), it is natural to consider the corresponding 
combinatorial Yamabe problem, i.e. does there exist any circle packing metric r with 
constant curvature Ri ? Inspired by [HI [53 [2a E3], we study the combinatorial Yamabe 
problem by the combinatorial curvature flows introduced in the following. 


3 The 2-dimensional combinatorial Ricci flow 


3.1 Definition of 2-diniensional combinatorial Ricci flow 


Ricci flow was introduced by Hamilton |29| to study the low dimension topology. For a 
closed Riemannian manifold the Ricci flow is defined as 


m 


9ij — 2i?j 


(3.1) 


with normalization 


3 2 

~ ‘2‘Rij-, 


(3.2) 


where r is the average of the scalar curvature. Ricci flow is a powerful tool and has lots 
of applications. Specially, it has been used to prove the Poincare conjecture @01 HD [12]. 
For a closed surface {M‘^,gij), the Ricci flow equation (3.1) is reduced to 




^Qij 


(3.3) 
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with the normalization 


(3.4) 


m 


gij = {r - R)g, 


IJ-, 


where R is the scalar curvature. It is proved [301 [I3] that, for any closed surface with 


any initial Riemannian metric, the solution of the normalized Ricci flow (3.4) exists for all 


time and converges to a constant curvature metric conformal to the initial metric as time 
goes to inhnity. It is further pointed out by Chen, Lu and Tian [12] that the Ricci flow can 
be used to give a new proof of the famous uniformization theorem of Riemann surfaces. 
The discrete version of surface Ricci flow was first introduce by Chow and Luo in [TT], in 
which they gave another proof of Thurston’s existence of circle packing theorem. In this 
subsection, we will introduce a different combinatorial Ricci flow on surfaces to study the 
constant curvature problem of Ri. 

Definition 3.1. For a weighted triangulated surface (M, T, 4>) with circle packing metric 
r, the combinatorial Ricci flow is defined as 


dgi _ „ 

dt ~ 


(3.5) 


where gi = r, 


(3.5) 


Following Hamilton’s approach, we introduce the following normalization of the flow 

dgi 


dt 


— {Rav Ri)9ii 


(3.6) 


where Rav is the average curvature of Ri defined by (2.3). It is easy to check that the 


tota l measure g{M) = ||r|p of M with respect to g, is invariant along the normalized flow 


(3.6), from which we know that the average curvature Rav = is invariant along 

(XB). As our goal is to study the existence of constant combinatorial curvature metric. 


we will focus on the properties of the flow (3.6) in the following. We assume r(0) G ^ 
and then r{t) G along the flow in the following. 


The flows (3.5) and (3.6) differ only by a change of scale in space and a change of 


parametrization in time. Let t,r,R denote the variables for the unnormalized flow (3.5), 


and t,r,R for the normalized flow (3.6). Suppose r{t), t G [0, T), is a solution of (3.5). 
Set r{i) = (/?(t)^/^r(t), where ip{t) = ||r||“^ and t = ip{T)dT. Then we have 

||f|p = 1, Ri = ip{t)~^Ri, Rav = (p{t)~^Rav = 2 -kx{M). 


This gives 

dgi 

di 
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where, in the last step, we use 




27rx{M)ip^. 


Conversely, if r{i),i € [0,T), is a solution of (3.6), set r{t) = where 


(p{i) = e ^ — I ^{T)dT. 


Then it is easy to check that ^ = —Rigi. 

For simplicity, we will set Ui = In gi in the following, 
could be written as 


Then the flows (3.5) and (3.6) 


ft = —R 


(3.7) 


and 

ft = Rayl - R (3.8) 

respectively, where u = (tti, • • • , un)'^, 1 = (1, • • • , 1)"^ and R = (i?i, • • • , Rn)'^- 


3.2 Discrete Laplace operator 

Since g is an analogy of the area element, we can define an inner product (•, •) on (M, T, d>) 
with circle packing metric r by 

N 

{f^h) = Y,fihr‘f = h^T.f (3.9) 

i=l 

for any real functions f,h G C'(C), where 

S = diag{rl,--- ,r^}. 

A combinatorial operator S : C{V) —)■ Ciy) is said to be self-adjoint if 

{Sf,h) = {f,Sh) 


for any /, /t G C(V}. 

The classical discrete Laplace operator m “A” is often written in the following form 

~ /O) 

i^j 

where the weight Wjj can be arbitrarily selected for different purpose. Bennett Chow and 
Feng Luo [H] first gave a special weight which comes from the dual structure of circle 
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patterns, while Glickenstein |24j gave a similar weight in three dimension. The curvature 
Ki is determined by the circle packing metric r. Consider the curvature as a function of 
u, where Ui = In gi = 2 In r* are coordinate transformations, then Bennett Chow and Feng 
Luo’s discrete Laplacian A.cl can be interpreted [20] as the Jacobian of curvature map, 
i.e. Acl = —2L, where 


L = 


- ,Kn) 


(3.10) 


d{ui, ■ ■ ■ ,un) ' 

It is noticeable that symbol L above is different from that in m by a factor 2, which 
comes from the coordinate transformations. 

Using the discrete measure /r, = rf, we can give a new definition of discrete Laplace 
operator, which is slightly different from that in [HEol. 


Definition 3.2. For a weighted triangulated surface (M, T, d>) with circle packing metric 
r : U —)• (0, +oo), the discrete Laplace operator A : C{V) —)■ CiV) is defined as 


Ah 




(3.11) 


for / e C{V). 

The following property of L will be used frequently in the rest of the paper. 


Lemma 3.3. ([Hj, Proposition 3.9.) L is positive semi-definite with rank N-1 and kernel 
{tl\t G IR}. Moreover, > 0, < 0 for f ~ j and = 0 for others. 

Using this lemma, we can write discrete Laplace operator as a matrix, that is. 


A = 


with 

A/ = -S-^L/. 

Note that the term — > 0 for i ~ j in the definition of the operator A, thus it could be 

taken as a weight on the edges. We sometimes denote this weight as Wij in the following 
if there is no confusion. This implies that A is a standard Laplacian operator defined on 
the weighted triangulated surface (M, 7~, <h) with circle packing metric r and measure /r. 
And it is easy to get 

{Af,g) = 5'^S(-E-1L/) = -g^Lf = {f,Ag). 

Hence the Laplacian operator A is self-adjoint with respect to (•, •). Moreover, if we define 

* j^i 
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then we can derive 


N 


(|V/|2,1) = ^ iV/LVf = Y,i-Lv){fj - fi? = fLf = -if, A/). 


(3.12) 


2 = 1 


J~2 


Using the discrete measure /r introduced in (2.1), (3.12) can be written in an integral form 


/ iV/pd^ = - / fAfd^i. 
Im Jm 


When the circle packing metric evolves along the flow (3.6), so is the curvature Ri. 


The evolution of Ri is very simple and it has almost the same form as the evolution of 
scalar curvature along the Ricci flow on surfaces derived by Hamilton in |3r)j . 


Lemma 3.4. Along the normalized combinatorial Ricci flow (3.6), the curvature Ri sat¬ 
isfies the evolution equation 


^ = AR, + R,{Ri - Ra 


Proof. As Ui = Inr?, we have 


I' 

dRi 

dt 


^ = Then we have 


E dRi duj 
duj dt 

Y{ — —-—r r -6- )(R - R ) 

2 o„,. A ' t')K^av 

. r- uuj r- 

4 -Rj)- R^{Rav - Ri) 

j 


~ r'^'^ f)n, T RiiRi Ra 


4 y.i-^)(Rj - Ri) + Ri(Ri - R-v)- 

^ ou 


r. 




'j 


In the last two steps. Lemma 3.3 is used. 


(3.13) 


Thus the evolution equation (3.13) of Ri is a reaction-diffusion equation. In fact. 


suppose H is a subset of U, set Vijf = fj — fi for / G C{V) and i ~ j, and denote 

dQ = {{ij] G E\i eUJ G O'"}, 

then it is easy to check that 


Then we have similar explanation as the smooth Ricci flow for equation (3.6). 
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3.3 Maximum principle 


As the evolution equation (3.13) is a heat-type equation, we have the following discrete 
maximum principle for such equations. Such maximum principle is almost standard, for 
completeness, we write it down here. 

Theorem 3.5. (Maximum Principle) Let f ‘-V x [0, T) — )• R be a function such that 

A/,+ $,(/,), yii,t)GVx[0,T) 

where <l>j : R —?• R is a local Lipschitz function. Suppose there exists Ci G R such that 
/i(0) > Cl for all i £ V. Let y? be the solution to the associated ODE 

/ S = 

I 9(0) = Cl. 


then 


fiit) > fit) 


for all (i, t) £ V X [0, T) such that fit) exists. 

Similarly, suppose / : E x [0, T) —)■ R be a function such that 

A/, + cl>,(/,), V(i,t)GPx[0 ,r). 

Suppose there exists (72 G R such that /j(0) < C 2 for all i £ V. Let V’ be the solution to 
the associated ODE 

\ >^(0) = C 2 , 

then 

fiit) < tpit) 

for all (i, t) £ V X [0, T) such that ^fit) exists. 


Remark 3. In fact, Theorem 3.5 is valid for general Laplacian operators defined as 

= '^aijit)ifj - fi), 




where atj > 0, but not required to satisfy the symmetry condition aij = a 




The proof of Theorem 3.5 is almost the same as that in m, we give a proof of it here 
just for completeness. 
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Proposition 3.6. Let / : P x [0, T) — )■ be a solution to the heat equation 

df 


dt 


= A/, 


i.e. 


dfi 

dt 


^fi = ^aij{t){fj - fi). 




If there are constants Ci < (72 G R such that (7i < /j(0) < C2 for all i G V, then 

Ci<m<C2, y{i,t) gv x[o,T). 

This proposition follows from the following more general theorem. 

Theorem 3.7. Suppose / ; P x [0, T) —)■ be a function with /j(0) > a for some 

dfi 


constant a G R, and / satisfies 


dt 


>Mi 


at any (i,t) gVx [0, T) such that fi{t) < a. Then 

fi{t) > a 

for all (z,t) G fo X [0,T). 

Proof. Note that H -.V x [0, T) —)■ R is a (7^ function and (i, to) is the vertex and 
time where Hi{t) attains its minimum among all vertices and earlier times, namely 


iLi(to) = min iT(t), 
(i,t)eVx[o,to] 


then 


dHi 

dt 


(to) <0,A77i(to) >0. 


Consider a function H defined by Hi{t) = (/i(t) — a) + et + e, where e is any positive 
constant. Then iLi(O) > e > 0, and 

dHi _ dfi ^ A . , _ A rr I 

~ “TT + ^ — ^fi + £ — ALtj + e 
dt dt 

for (z,t) G P X [0,T) satisfies fi{t) < a. Thus we just need to prove that Hi{t) > 0 for 
(i,t) GPx [0,T). 

Suppose that (f,to) is the first point in P x [0,r) such that iLj(to) = 0. Then 


Then we have 


/i(to) = Hi{to) + a - eto - e = a - eti) - e < a. 


OH 

0 > > AiL*(to) + e > e > 0, 


which is a contradiction. And the theorem follows from e —)• 0. 
Now suppose (3 :V x [0, T) —)■ R be a function. 
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Proposition 3.8. Let f '■ V X [0, T) —)• R be a function such that 

A/, + A/,. 

And for each r G [0, T), 3Cr < oo such that Pi{t) < Ct for (i, t) G P x [0, r]. If /i > 0 for 
all i €V, then fi{t) > 0 for all {i,t) G P x [0, T). 

Proof. Given r G (0,r), define 

Ji{t) = 


then 


dJj 

iH 


> AJ, + (A - Cr)Ji. 


Since /3i(t) < Cr for all {i,t) £V x [0,r], then we have 

dJi 


dt 


for {i,t) G P X [0,r] such that Ji{t) < 0. Then by the Theorem 3.7, we have Ji{t) > 0, 
and then fi{t) > 0 for any {i,t) G P x [0,T). ■ 

Now we give the proof of Theorem |3.5[ 


Proof of Theorem 3.5t We take the lower bound to prove, and proof for the upper 
bound is similar. By assumption, we have 


-<p)> -‘f) + ^iifi) - 

The assumptions on the initial data imply that 

(/*-V^)(0) >0,ViGP. 

Let r G (0, T), then 3Cr < oo such that \fi{t)\ < Cr for all (i, t) G Px [0, r] and |(^(t)| < Cr 
for t G [0,r]. Since F is locally Lipschitz, then we have 

\Mfi)-^i{^)\<Lr\fi-ip\,yieV 

for some positive constant Lr > 0. Then 


^{fi -<p)> -‘f)- LrSgn{fi - (f) ■ {fi - G P x [0,r]. 


Applying Proposition |3.8[ we get 

fi{t) - Lp{t) > 0,V(i,t) G P X [0,r]. 

And the theorem follows from r G [0, T) is arbitrary. ■ 

Applying Theorem 3.5 to (3.13) with vanishing initial value, we can easily get the 
following corollary. 
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Corollary 3.9. If Ri{0) > 0 (-Rj(O) < 0) for all i £ V, then Ri{t) > 0 {Ri{t) < 0) for all 
i £ V as long as the flow exists, i.e. the positive and negative curvatures are preserved 


along the normalized flow (3.6). 


Set Rmax{t) = maxjgyi?i(t) and Rmin{t) = uiini^v Ri{t)- Applying Theorem 3.5 


to 


(3.13) with general initial data, we get the following lower bounds for Ri{t). For simplicity, 


We just state it in the form parallel to that given in [El- 

Lemma 3.10. (Lower Bound) Let ri{t) be a solution to the normalized combinatorial 


Ricci flow (3.13) on a closed triangulated surface (M, T, <I>). 
(1) If x(M) < 0, then 


Ri — Rav > 


Ra 


_ ^ay _^ pRayi 

^min ( 0 ) 


1 — (1 — 


Rav ^ (Rmin(O) Rav)(^ 


Rav^ 


(2) If x(M) = 0, then 


Ri > 


^min ( 0 ) 


1 

> 


1 - Rmini0)t t 

(3) If x(M) > 0 and Rmin{0) < 0, then 

Rav 


R > 


1 -( 1 - 


\pRav^ 


^ -Rmm( 0 )c 


Ravi 


Notice that in each case, the right hand side of the lower bound estimate tends to 0 as 
t —)■ oo. However, for the upper bound, the situation is not so good. The solution of the 
ODE 

/ S = s(s — Rav) 

I s{0) = So 


corresponding to (3.13) is 


s{t) = < 


0 , 

so 
l-soi ’ 


Ra 


1 — (1 — ) e-Rat) t 


So = 0 

so / 0, x(M) = 0 
, -So 7 ^ 0, x(M) ^ 0 


If So > max{i?aD,0}, there is T < oo given by 

T = 


si. xiM) = 0 
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such that 


lim s{t) = + 00 . 
t^T- 


The implies that, in the case of Rmax{0) > max{i?a^,0}, Theorem 3.5 will not give us 
good upper bounds for the curvature along the combinatorial Ricci flow (3.13). However, 
in the case of RmaxiO) < 0, we have good upper bounds. 


Lemma 3.11. If i?i(0) < 0 for all i G H, then we have 

Ri{t) - Ray <Ray(l- ' 

V -^maxiy) J 

These results are almost parallel to that of the smooth Ricci flow on surfaces. 


3.4 Convergence of the 2-dimensional combinatorial Ricci flow 


Combining Lemma (3.10) with (3.11) gives 


Theorem 3.12. Suppose (M, 7~, 4>) is a weighted triangulated surface with initial circle 
packing metric r(0) satisfying i?j(0) < 0 for all i G V, then there exists a negative con¬ 
stant curvature metric r* on (M, T, d>). Furthermore, the solution r(t) to the normalized 
combinatorial Ricci flow ( |3.6[ ) exists for all time and converges exponentially fast to r* as 
t —)• -|-oo. 


Theorem 3.12, derived by combinatorial maximum principle, implies more than it 


seems. It claims that, if there is a metric with all curvatures negative, then there always 
exists a negative constant curvature metric, and vice versa. Using this fact, we will give 
a combinatorial and topological condition which is equivalent to the existence of negative 
constant curvature metric in subsection 3.6 For the convergence of combinatorial Ricci 
flow, the negative initial curvature condition in Theorem 3.12 is a little restrictive. It is 
natural to consider Hamilton’s approach to generalize this result. However, we found that 
there is some technical difficulties to go further with Hamilton’s approach in this case. By 
studying critical points of (3.6), which is an ODE system, we can get the following local 
convergence result. 


Theorem 3.13. Suppose r* is a constant R-curvature metric on a weighted triangulated 
surface (M, T, <I>). If the first positive eigenvalue Ai of —A at r* satisfies 


ai(-a)>r:,= 


2vry(M) 


(3.14) 


and ||r(0) — is small enough, then the solution to the normalized combinatorial Ricci 


flow (3.6) exists for t G [0, -|-oo) and converges exponentially fast to r*. 
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Proof. We can rewrite the normalized combinatorial Ricci flow (3.6) as 


dt - 2^^^^ 

Set rj(r) = \{Rav — Ri)ri, then the Jacobian matrix of r(r) is given by 


d /I 


{DrT{r))ij = — ( -{Rav - Ri)ri ) = -(Rai; “ Ri)Sij + {Ri^ij “ R 


1 


nr 




1 dKj 

TiVj dUj 


Denote A = S aLS 2 ^ then A ~ S 2 AS 2 = —A, hence Ai(—A) = Ai(A). At the 
constant curvatnre metric point r*, we have 

r.T 


DrFL* = Rav I - 


rr 


— S 2 AS 2 = R 


I - 




-A. 


Select an orthonormal matrix P such that 

AP = diag{0, Xi{A),-■ ■ ,AAr_i(A)}. 

Suppose P = (eo, ei, • • • , cat-i), where e* is the {i + l)-column of P. Then Aeo = 0 and 
Aci = AjCj, 1 < f < Ai — 1, which implies cq = ''’/||?'|| and r T e^, 1 < f < A^ — 1. Hence 
{lN-j^)eo = 0 and {lN-j^)ei = Ci, 1 < i < N -1, which implies P'^{lN-j^)P = 
diag{0, 1, • • • ,1}. Therefore, 

= P ■ diag{0,Rav - Ai(A), • • • ,Rav - AAr_i(A)} • P'^. 

IfAi(A) > then DrT\^, is negative semi-definite with kernel {cr|c G M} and 

rank (Drr|^*) = Ai — 1. Note that, along the flow (3.6), ||r|p is invariant. Thus the kernel 
is transversal to the flow. This implies that DrT\r* is negative definite on and r* 

is a local attractor of the normalized combinatorial Ricci flow (3.6). Then the conclusion 
follows from the Lyapunov Stability Theorem! |43j. Chapter 5). ■ 


If x(M) < 0, we always have Ai(—A) > 0 > R*^ = thus we get 

Corollary 3.14. Suppose there is a nonpositive constant curvature metric r* on a 
weighted triangulated surface with x(M) < 0. If ||r(0) — r*|p is small enough, 

then the solution to the normalized combinatorial Ricci flow ( |3.6| ) exists for t G [0, -|-oo) 
and converges exponentially fast to r*. 


Before proving the global results, we first check that the convergence of the flow (3.6) 
ensures the existence of the constant curvature metric. More generally, we have the fol¬ 
lowing result. 
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Lemma 3.15. Suppose that the solution to the flow (|3.6|) lies in a compact region in 
IR^g, then there exists a constant curvature metric 
of metrics which converge to r* along this flow. 


Moreover, there exists a sequence 


Proof. Consider the Ricci potential 

P{u) = / '^{Ki - Ravrf)dui, 

Juo 

where Ui = Inr? and uq G is an arbitrary point. If we set ipi = Ki — by 

direct calculations, we have 

d(pi _ dipj 
duj dui 

As the domain of u is simply connected, this implies that the integration is path indepen¬ 
dent and then the Ricci potential is well defined. Suppose t G [0, T) and T is the maximal 
existing time of u{t). Then {rt(t)} CC implies T = -|-cx), otherwise u{t) will run out 
of the compact region. {ri(t)} CC also implies that F{u{t)) is bounded. Moreover, 

i-^FW)) = (V„F)^ = - E(-R. - < 0. 

\ * / j 

Hence limi_^_|_oo F{u{t)) exists. Then by the mean value theorem, there exists a sequence 
^ in,n + 1) such that 

F{u(n + 1)) - F{u{n)) = (F(n(t)))'|i=^„ = - - RavfrfiU) ^ 0 

i 

as n —> -|-oo. Since {r(f)} CC R> 0 ) further select a subsequence tk = ink such 

that r{tk) r* as k ^ -|-oo. Combining Ylir‘f{Ri — Rav)‘^{ink) 0; we know that r* 
determines a constant R-curvature. ■ 


Corollary 3.16. Suppose that the solution to flow (3.6) exists for all time and converges 


to r(-|-oo), then r(-|-oo) is a constant curvature metric. 


Remark 4. For the Ricci potential functional F used in the proof of Lemma 3.15 the 
primitive form of this type of functional was first constructed by Colin de Verdiere in m 
and then further studied by Chow and Luo in HI- Our definition is a modification of 
their constructions. 


For the Ricci potential F introduced in the proof of Lemma 3.15 we further have the 
following property. 
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Lemma 3.17. Given a weighted triangulated surface (M, T, <h). Assume there is a con¬ 
stant curvature circle packing metric r* on Define the Ricci potential as 


N 


^(«)= / - Ra,rf)d 


Ui. 


(3.15) 


i=l 


Denote 


— {u £ Ill'll Ui = a}, a G R. 


If Ai(—A) > Rav for all r G 

(1) HesSuF is positive semi-definite with rank N — 1 and kernel {tl\t G R}. 

(2) Restricted to F\>^^ is strictly convex and proper. has a unique zero point 

which is also the unique minimum point. Moreover, lim F(u) = -|-oo. 

u£‘Wa,U^OO 

Proof. By direct calculations, we have 


Hess,,F = L — Rnv'^^ ( I — tt— r? 


Sa = S 2 ( A - R„ 


I - 


rr 


1 

Sa. 


(3.16) 


By the same analysis as that of Theorem 3.13, we know that, if Ai(—A) > Rav, HesSuF 
is positive semi-definite with rank{HesSuF) = N — 1 and Ker{HesSuF) = {tl|t G R}. 

For the second part of the proof, we follow that of Colin de Verdiere m and Chow and 
Luo [IT]. It’s easy to see F{u) = F{u -|- tl) for any t G R and u G R'^. As F is invariant 
along the direction 1, we just need to prove it on = {u £ R'^l — 0 }- ^ rigorous 


proof is formulated in Theorem B.2 of 
point and minimum point of F\i^^. 


In addition, u* + 


a-l^ u* 
N 


1 is the unique zero 


Theorem 3.18. Suppose is a weighted triangulated surface and Ai(—A) > 

Rav for all r £ R>o- Then the solution to the flow (3.6) converges if and only if there 
exists a constant curvature metric r*. Furthermore, if the solution converges, it converges 
exponentially fast to the metric of constant curvature. 


Proof. The necessary part can be seen from Corollary 3.16 For the sufficient part, 
without loss of generality, we may assume r* £ Consider the flow (3.6) with initial 

metric r(0) G S^“^nR>o- Then r{t) £ S-^^^nR^g and u{t) £ Tj* = {u G R^| Yli = 1} 
for all t. Let 11 be the orthogonal projection from to the plane = {u £ R'^l Yli = 
0}. Then F{u) = T(n(u)) and n(u) —)■ oo as u —)■ oo. In fact, for a sequence 
n(u("’)) is unbounded if and only if |n(u("’))j — n(u("'))j| = is unbounded. 

lim F(u) = - 1 - 00 , (3.17) 

u^oo,u&'r* 


Then Lemma 3.17 implies 
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which means that F{u) is still proper on Fi . Since F(u{t)y < 0, u{t) lies in a compact re¬ 
gion in yy. By Lemma pl.15 the solution exists for t G [0, -|-oo) and there exists a sequence 
of metrics r(tfc) which converge to r* as tk t +oo. Hence F{u{+oo)) = lim F{u{tk)) = 

_ fc^+OO 

F{u*) = 0. While Lemma 3.17 says that u* is the unique zero and minimum point of F, 
thus we get u{t) converges to u* as t —)• -foo. The exponentially convergent rate comes 
from Theorem 13.131 ■ 

Corollary 3.19. Suppose is a weighted triangulated surface with y(M) < 0. 

Then the solution to the flow ( |3.6[ ) converges if and only if there exists a constant R- 
curvature metric r*. Furthermore, if the solution converges, it converges exponentially 
fast to the metric of constant curvature. 


We expect Corollary 3.19 to be still valid for surfaces with x(M) > 0. However, things 
are not so satisfactory for general triangulations. In fact, things may become very different 
and complicated. Assume there exists a constant curvature metric r*, we may compare 
the Ricci potential (3.15) with 


Giu)= / Y.^Ki-Kl)du^, 

i=l 


where K* = then 


N 


F{u) - G{u) = r '^{K* - Ravrf)dui 

- u*) 


2 = 1 


= 27rx(M) 


= 27rx(M). 


11 2 


e^^dui + • • • -I- e^^duN 

giii e“jv 




112 


- 27 rx(M) ln(e“i +--- + e“^) 


Restricted to the hypersurface F* = {u ^ 1R'^| term becomes 

zero, hence 


F{u) = G{u) + 27rx(M)- 


1^^ 112 


Vrt G y* 


By similar arguments in the proof of Theorem |3.18 , we further have 

lim G{u) = + 00 . 

u^oo^uG'f* 

One may expect that the growth behavior of F is similar to that of G. Unfortunately, 
there is an obstruction which makes things very complicated. It’s easy to see 




lim 

u^oo^uG'f* 


11 2 


= —oo. 
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The growth behavior of F is uncertain unless one can compare the growth rates of these 
two terms concretely. It seems that we can not expect that G{u) succeeds the last term 
due to the following example. 


Example 2. Given a topological sphere, triangulate it into four faces of a single tetra¬ 
hedron. Fix weights = 0. Denote r* = 1 , then r* is a constant curvature metric with 
Ri = vr. By direct calculation, we get HesSuF\u* = — f )(4/4 — 11^), which is negative 

semi-definite with three negative eigenvalues. Up to scaling, HesSuF is negative definite 
at r*. Thus the fixed point r* = 1 of the ODE system = ^{Rav — Ri)F (i.e. the flow 
(3.6)) is a source. For this particular weighted triangulation, the flow (3.6) can never con¬ 
verge to the constant curvature metric r* for any initial value r(0), as t —)• -|-oo. However, 
—oo, the solution of (3.6) converges to r* if r(0) G §'^“^(2) is close enough to 


when t 
r . 


From the example we know that it is impossible to get an analogous version of Corollary 


3.19 when x(M) > 0. However, we can get the following long time existence of (3.6) by 


maximum principle. 

Theorem 3.20. Given a weighted triangulated surface (M, T, d>) with xi^) > 0 

(tt - $(e)) + 27ry(F7) < 0, V/: 0 g / C y. ( 3 .I 8 ) 


E 

{e,v)&Lk{I) 


Suppose ro is a circle packing metric with Ri > 0 for all i G U. Then the solution of (3.6) 
with initial metric tq exists for all time and lies in a compact region in H^g. Furthermore, 
there exists tn t + 00 , such that r(t„) converges to a constant curvature metric r*. 

we just need to prove that {r(t)} CC H^g. We prove it by 


3.15 


Proof. By Lemma 
contradiction. Suppose 

Mt)} n 9(]R^g n ( 

then there exists a sequence tn such that 

r{tn) —>• r* G (9(R>g PI 

where r* = (0, • • • ,0, ’’J/i+i 
sition 4.1 of |14j . we have 


lim y Ki{r{tn)) = - 
n^+oo 

i£l {e,v)GLk{I) 


r^) for some nonempty proper subset / of U. By Propo- 
^ (tt - 4>(e))-h 27rx(T>) < 0. (3.19) 


However, by Corollary 3.9, Ri > 0 is preserved along the flow (3.6). Thus Ki > 0 for all 
i G U, which implies lim„_,._|_oo This contradicts (3.19). ■ 

Remark 5. It is easy to see that the triangulation of sphere in Example does not satisfy 
the condition (3.18). 
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3.5 The prescribing curvature problem 

Using modified combinatorial Ricci flow, we can consider the prescribing curvature problem 
on a weighted triangulated surface (M, T, d>). 


Definition 3.21. Suppose (M, 7~, <h) is a weighted triangulated surface with circle packing 
metric r, R € C(y) is a function defined on M. The modihed combinatorial Ricci flow 
with respect to R is defined to be 


where gi = rf as before. 


dgi 

dt 


{Ri 


(3.20) 


Note that the total measure g-{M) = ||r|p may vary along the modihed combinatorial 
Ricci how (3.20), which is different from that of (3.6). 

R is called admissible if there is a circle packing metric r with curvature R. For a 
given function R G C(V), we can introduce the following modihed Ricci potential 


Fiu) 



(3.21) 


It is easy to check that the modihed Ricci potential F is well-dehned. Furthermore, by 
direct calculations, we have 


/ Ri 


HesSniF = L — 


1 

S2. 


Rn 


The following lemma is useful. 

Lemma 3.22. ([ 13 ]) Suppose C IR'^ is convex, the function /i : —>■ K, is strictly 

convex, then the map Vh : If —)• IR-^ is injective. 


Theorem 3.23. Suppose (M, T, $) is a weighted triangulated surface and R G C{V) is 
a function dehned on M. 


( 1 ) 

( 2 ) 


If the solution to the modihed how (3.20) converges, then R is admissible. 


If Rj < 0 for all i, but not identically zero, and R is admissible by a metric r. Then 
r is the unique metric in R^q such that its curvature is R. Moreover, the solution 
to the modihed how (3.20) converges exponentially fast to r. 
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Proof. The first part is obviously, and R is admissible by metric r(+oo). For the 
second part, notice that 

HesSuF = L — T,2diag{Ri, ■ ■ ■ ,i?Ar}S2, 

it is easy to check that, if i?, < 0 for i = 1, • • • ,N and not identically zero, HesSuF is 
positive definite. By Lemma 3.22, = (iFi — Rirl, ■ ■ ■ ,K]\[ — RNr^^)"^ is an injective 

map from u G to R^. Hence f is the unique zero point of VuF. This fact implies that 
r is the unique metric in R^q such that it’s curvature is R. By Lemma B.l in |20j . we know 
that F is proper and lim F{u) = +oo. Furthermore, ^F{u{t)) = — {Ki — Rirf^ < 


0 implies that the solution to (3.20) lies in a compact region. The rest of the proof is the 


same as that of Theorem 3.18, so we omit it here. 


Remark 6. Theorem 3.23 could be taken as a discrete version of a result obtained by 
Kazdan and Warner in m- 


Remark 7. The second part of Theorem 


3.23 


i, the corresponding prescribing curvature problem is already solved in Corollary 3.19 


implies that x(M) <0. If = 0 for all 

In 


this case, the metric r is not unique. However, it is unique up to scaling. This is slightly 
different from Theorem 13.231 


3.6 Existence and uniqueness of constant curvature metric 

We first condider the uniqueness of constant i?-curvature metric. 

Theorem 3.24. Suppose (M, 7~, <h) is a weighted triangulated surface, c* G R is a 
constant. 


(1) If c* < 0, there exists at most one circle packing metric with curvature R = c*l. 

(2) If c* = 0, there exists at most one metric with curvature R = c*\ up to scaling. 


Proof. The first part is just a corollary of the second part of Theorem 3.23 The 
second part is already proved by Thurston 


If x(M) < 0, Theorem 3.24 implies that the constant curvature metric (if exists) is 
unique or unique up to scaling. If x{FI) > 0, there may be several different constant 
curvature metrics on and we have the following example. 


Example 3. Consider the weighted triangulation of the sphere in Example with tan¬ 
gential circle packing metric. Let ri = 1, r* = x > 0, 2 < f < 4 and then x = 
where 9 is the inner angle of AV 1 V 2 V 3 at the vertex U 3 . 
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By direct calculations, Ki = 69 — tt, Ki = — 26, 2 < i < A. Constant i?-curvature 

implies that 


|7r-20 


60 — TT 


Using MATLAB, we can solve the equation and get two solutions x = 1 and x 
5.9487. Thus we find a 
(1,5.9487,5.9487,5.9487). 


5.9487. Thus we find at least two constant curvature metrics r\ = \ and 


In the following of this subsection, we consider the existence of constant 72-curvature 
metric. For the classical discrete Gauss curvature Ki, Theorem |l.l| states that the existence 
of constant 77-curvature metric is equivalent to 


27rx(M)|^ > - ^ (tt - 4>(e))-7 27rx(F/) 

' ' {e,v}£Lk(I) 


(3.22) 


for any nonempty proper subset 7 of U. Theorem 1.2 further states that the space of all 
admissible classical discrete Gauss curvature K is 


^ 4 


j^gb n ( n 


where A^gb and are defined by (1.6) and <113 respectively. If r* determines a constant 
72-curvature, then the 77-curvature 77* = (Kf,--- is admissible, where 77* = 

Hence 77* G 'S'l, i.e. 




2 vry:(M)^p^^^ > - ^ {tt - <^{e)) + 2TrxiFi) 


(3.23) 


{e,v)&Lk{l) 


for any nonempty proper subset 7 of V. (3.23) is necessary for the existence of constant 


72-curvature metric, which involves the combinatorial, topological and metric structure of 
the triangulated surface. In the case of x(M) < 0, the combinatorial-topological-metric 


condition (3.23) is also sufficient. In fact, we have the following result. 



Figure 2: tetrahedron triangulation of sphere 
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Theorem 3.25. Suppose {M,'T,^) is a weighted triangulated surface with x{^) ^ 0- 
Then there exists a constant i?-curvature metric if and only if there exists a circle packing 


metric r* such that, for any nonempty proper subset I of vertices V, (3.23) is valid. 


Proof. In the case of xi^) = 0, a zero i?-curvature metric is exactly a zero in¬ 
curvature metric. The conclusion is contained in Theorem ll.il 

In the case of x{M) < 0, We just need to prove the sufficient part. Set K* = 


27rx(M) 




By Theorem 


1.2 


(3.23) implies that K* = ,in^)^ is an ad¬ 


missible in-curvature. Suppose K{rQ) = K*. Then rg is a circle packing metric with 


Ri{fo) < 0 for all i (z V. By Theorem |3.12[ the combinatorial Ricci flow (3.6) with tq as 
the initial metric converges to a negative constant ii-curvature metric. Therefore, there 
exists a constant i?-curvature metric on (M, T, <h). ■ 

The key point of the proof of Theorem |5.7| is to find a circle packing metric with all 


ii-curvature negative and then use the combinatorial Ricci flow (3.6) to evolve it to a 


constant ii-curvature metric. Therefore, we can generalize Theorem 5.7 to the following 
form. 


Theorem 3.26. Suppose (M, T, <h) is a weighted triangulated surface. 

(1) There exists a negative constant curvature metric if and only if n 1 R:^q ^ cj)- 

(2) There exists a zero curvature metric if and only if 0 G ‘3^. 

(3) If there exists a positive constant curvature metric, then 3^ n R^q ^ cj). 


We include the case of x{^) > 0 in Theorem 3.26 In this case, if there exists a 


constant ii-curvature metric, then there exists a circle packing metric r* satisfying (3.23), 
which implies 3^ n R^q / 4 >. 


Similar to Thurston’s condition (1.4), conditions in Theorem (3.26) also show that the 


combinatorial structure of the triangulation and the topology of surface, which have no 
relation to the circle packing metrics, contain some information of R-curvature. 

For the surfaces with x{^) > 0, the property of constant R-curvature metric seems 
complicated and confusing as exhibited in Example and Example We do not know 
whether the condition 3/^ n R^q 7 ^ </) is sufficient for the existence of constant R-curvature 
metric. Even through, we can get some sufficient conditions for the existence of constant 
R-curvature metric by the discrete maximum principle. 


Theorem 3.27. Suppose (M, T, 4*) is a weighted triangulated surface. If there is a circle 
packing metric r satisfying Rj > 0 for all z G R, and 

^ (tt - 4>(e))27rx(R/) < 0 

{e,v)^Lk{I) 
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for any nonempty subset I oi V, then there exists a non-negative constant i?-curvature 
metric on (M, T, ‘h). 


Proof. By the maximum principle, i.e. Theorem 3.5 > 0 is preserved along 


the normalized combinatorial Ricci flow (3.6). Using similar arguments in the proof of 


Theorem 3.20 we get the proof. 


4 Combinatorial Calabi flow on surfaces 


Given a compact complex manifold admitting at least one Kahler metric, to find the 
extreme metric which minimizes the norm of the curvature tensor in a given principal 
cohomology class, Calabi [7] introduced the Calabi flow, which could be written as 




(4.1) 


on a Riemannian surface. Chrusciel |16j proved that the Calabi flow (4.1) exists for all 
time and converges to a constant Gaussian curvature metric on closed surfaces using the 
Bondi mass estimate, assuming the existence of the constant Gaussian curvature in the 
background which is ensured by the the uniformization theorem. Chang m pointed 
out that Chrusciel’s results still hold for arbitrary initial metric, which implies the uni¬ 
formization theorem on closed surfaces with genus greater than one. Chen m gave a 
geometrical proof of the long-time existence and the convergence of the Calabi flow on 
closed surfaces. The first author [20] first introduced the notion of combinatorial Calabi 
flow in Euclidean background geometry and proved that the convergence of the flow is 
equivalent to the existence of constant classical discrete Gauss curvature. Then we [22] 
studied the combinatorial Calabi flow in hyperbolic background geometry. We also use 
the combinatorial Calabi flow to study some constant combinatorial curvature problem 
on 3-dimensional triangulated manifolds [23j. In this section, we study the R-curvature 
problem by combinatorial Calabi flow. 

Definition 4.1. For a weighted triangulated surface (M, T, <h) with circle packing metric 
r, the combinatorial Calabi flow is defined as 

dt 

where gi = rf and A is the Laplacian operator given by (3.11). 


(4.2) 


It is easy to check that the total measure fi{M) = ||r|p of M is invariant along the 

Interestingly, the combinatorial curvature Ri evolves 


combinatorial Calabi flow (4.2) 
according to 


^ = -Ri^Ri - A^Ri, 

dt 
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which has almost the same form as that of the scalar curvature along the smooth Calabi 
flow on surfaces. We can rewrite the combinatorial Calabi flow (4.2) as 


d'Ti _ _-__ p i 

M ~ 2ri ^ ^ ~ 


2 ri 


Set 




If there exists a constant curvature metric r*, then we have 


DjFi\r* — 


E 


1 dKi dKj 
rl duk duk 


— R„ 


dKj 
’ duj 


If we further assume that the condition (3.14), i.e. Ai(—A) > = ^|pr|^, is satisfied, 

then nr I,.* is negative semi-definite with rank N — 1 and kernel {t\\t G R}. From this 


fact we can get local convergence results similar to Theorem 3.13 Corollary 3.14 We can 
also get results similar to Lemma 3.15, Corollary 3.16, Theorem 3.18 and Corollary |3.19 
We just state the following main theorem for this section here. 

Theorem 4.2. Suppose (M, T, 4>) is a weighted triangulated surface with x(M) < 0, 


then the combinatorial Calabi flow (4.2) converges if and only if there exists a constant 
i?-curvature metric r*. 


Proof. If the solution of (4.2) converges, u(-|-oo) must be a critical point of this ODE 


system, i.e. An(-|-oo) = —S ^LR{+qo) = 0, which implies that i?(-|-oo) belongs to the 
kernel of L and hence r(-|-oo) is a constant curvature metric. 


Assume there exists a constant curvature metric r*. Write the flow (4.2) in the matrix 


form ii = AiZ = —S ^LR or = —LR, where = (rf, • • • , It is easy to see that 

j^F{u{t)) = -{K - Ravr^f^~^LR = -R^LR < 0. 

Using the properties of the Ricci potential, we can get the convergence result by similar 
arguments in the proof of Theorem |3.18 ■ 


Analogous to the combinatorial Ricci flow, we can also use the combinatorial Calabi 
flow to study the combinatorial prescribing curvature problem. In fact, we have the 
following result. 
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Theorem 4.3. Suppose (M, T, $) is a weighted triangulated surface, R G C{V) is a 
function defined on M with Ri < 0 for all i. Then R is admissible if and only if the 
solntion of the modihed combinatorial Calabi flow 

^=A{R- R),gi (4.3) 

exists for all time and converges. ■ 


The proof is jnst a combination of that of Theorem 3.23[ Remark and Theorem 4.2 
we omit it here. 


Following [20], we can also introduce a notion of energy to study the combinatorial 
cnrvatnre Ri. 


Definition 4.4. For a weighted triangulated surface (M, T, <h) with circle packing metric 
r, the combinatorial Calabi energy is dehned as 

N 

Cir) = Y,^l (4.4) 

i=l 


where (pi = Ki — 


27IXM „2 
llrIP i ' 


Consider the combinatorial Calabi energy C as a function of u, we have V^C = 2 A'^p>, 
where 


A = 


(9(ui, • • • , un ) 


dyi 

dyi 

dui 

duj\f 

dVN 

dVN 

dui 

du]si 


Note that A is in fact given by (3.16). Thus, if x(M) < 0, ^4 is symmetric, positive semi- 
definite with rank N — 1 and kernel {tl|t G R}. Using this fact, we can dehne a new flow, 
which is the gradient flow of C. We call it the modified combinatorial Calabi flow. 


Definition 4.5. For a weighted triangulated surface (M, T, 4*) with circle packing metric 
r, the modified combinatorial Calabi flow is defined as 


u = -^VuC, (4.5) 

or equivalently, 

ii = —A^ip, (4.6) 

where (p = {(fi,--- ,(Pn)'^- 
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Note that YliV’i — Oj so ™ oase of x{M) < 0, if the flow (4.5) converges, it 


converges to the constant curvature metric. Along the flow (4.5), we have 


(f = -AA^ip, C = -2\\A^p>\\^ < 0 


and 


^F{u{t)) = .u = ip'^.u = -(p^A^ip < 0 , 


which implies that F{u{t)) is decreasing along the flow (4.5). 


Following the arguments in the proof of Theorem 3.18 and 4.2 we can derive the 
following result. 

Theorem 4.6. Suppose (M, T, 4>) is a weighted triangulated surface with x{^) ^ 0) 
then the existence of constant curvature metric r* is equivalent to the convergence of the 


modified combinatorial Calabi flow (4.5) 


5 Combinatorial ct-curvature and combinatorial ct-flows 


In the previous sections, we investigated the properties of the new discrete Gauss curvature 
Ri = Ki/r‘f. Since the area of the disk packed at i is just vrr?, the denominator of 
the curvature rf (omitting the efficient vr) may also be considered as an “area element” 
attached to vertex i. We want to know whether there exists other types of “area element”. 
Suppose Ai is a general “area element”, which is an analogy of the volume element in the 
smooth case, the combinatorial Gauss curvature could be defined as Ri = Ki/Ai. The 
average curvature should be 27rx(M)/^Aj, which is an analogy of f Rdp/ f dp in the 
smooth case. If we expect the functional 


F{u) 



2t:x{M) 


Ai)dui 


to be well defined, the following formula 


A 

dui 


Efc 


d 

duj 



(5.1) 


should be satisfied for all i and j. It is a trivial observation that Ai = rf always satisfies 
(5.1) for all a G R. 

When x(m) ^ 0) have already seen in the previous sections that the combinatorial 
Ricci flow and Calabi flow are good enough to evolve circle packing metrics to a constant 
R-curvature metric. However, in the case of x{^) > Oj do not know how to evolve 
it. Maybe this is because that the growth behavior of F is uncertain (in integral level), 
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or equivalently, HesSuF is not positive semi-definite (in differential level), even for the 
simplest triangulation of (see Example . 

The above two reasons motivate us to consider a new “area element” Ai = rf. We 
find that, for x(M) > 0, if a < 0, combinatorial flow methods are good enough to evolve 
the curvature to a constant. For x(M) < 0, if a > 0, the properties of combinatorial flows 
are almost the same with previous sections. It is very interesting that m and [20] can be 
included into the case of a = 0. 


Definition 5.1. For a weighted triangulated surface (M, T, <h) with a circle packing metric 
r, the combinatorial a-curvature is defined as 


Ki 


(5.2) 


where a is a real number. 


For this type of curvature, we can also consider the corresponding constant curvature 
problem and prescribing curvature problem. As the methods are all the same as that we 
dealt with i?-curvature, we will give only the outline and skip the details in the following. 
The measure defined on M is now ^a{i) = ff and the average curvature is now 

ad^a 27rx(M) 

^ f da ^ ||rll« ’ 

Jm M ' Ma 

where ||r||Q, = Ikllo is defined to be N in the case of a = 0. In this 

section, we set u = In and define the Ricci potential as 

PU ^ 

Fa{u) = / '^{Ki - Ra,avrf)dUi, 

Juo 


where uq G is an arbitrary point. By direct calculations, it is easy to check that the 
Ricci potential Fa is well-defined. We further have VuFa = K — Ra,avT'°^ and 


Hess,, Fa = L — aR, 


a,av 


(rl 

V 



^\T 


I - 


r 2 . (r 2) 


/ - 
/ r? 


V 


/ 


where r“ = (rf, • • • ,rfr)'^, Aq = S 4 LS 4 and Lij = with Uj = Inrj. Note 


dK, 


q; , a , 


that the matrix I — has eigenvalues 1 {N — 1 times) and 0 (1 time) and kernel 
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{cr2 |c G R}. Following the arguments in the proof of Lemma 
positive eigenvalue of Aq, satisfies 

^ 0'Ra,avi 

HesSuFa is positive semi-definite with rank N — 1 and kernel {cl|c G R}. Especially, if 
ax(M) < 0, HesSuFa is positive semi-definite with rank N — 1 and kernel {cl|c G R}. 

We can also define the following modification of the combinatorial Ricci flow, which is 
called a-Ricci flow. 


3.17 


we have, if the hrst 


Definition 5.2. For a weighted triangulated surface (M, T, <h) with circle packing metric 
r, the a-Ricci flow is defined to be 


dn 

— = -Ra,iri. 

at 


The normalization of the a-Ricci flow is 


— {Ra,av 


(5.3) 


(5.4) 


and the normalized flow (5.4) are just Chow and 

■N 


Notice that, when a = 0, the flow 
Luo’s combinatorial Ricci flows mj. In this case, (or ^i) is invariant along the 

normalized flow (5.4). When a 7^ 0, ||r||Q, (or e““*) is invariant along the normalized 

icci fl 

1 / dK^' 


flow (5.4). Along the normalized a-Ricci flow, the a-curvature evolves according to 

N 


dR, 


a,I 


dt 


E 

i=i 


dui 


Fa,j T OiRa^i {^RoL^i - Ra 


(5.5) 


For a = 2, this is just the evolution equation (3.13) derived in Lemma 3.4 As the 
evolution equation (5.5) is still a heat-type equation, the discrete maximum principle, i.e. 
Theorem 3.5, could be applied to this equation. The evolution of curvature (5.5) under 
the normalized a-Ricci flow suggests us to define the a-Laplacian as 

N 


= 4 E 


i=i 


OKi 


duj ) r: 


1 


E 

j^i 


OKj 

dui 


ifj-h), 


(5.6) 


where / G C{V) and Ui = Inr*. Using this Laplacian, we can define the combinatorial 
a-Calabi flow as 

^ = {AaRa)iri. 


(5.7) 


Similarly, denote ifa,i = Ki — 
a-Calabi flow as 


27rx(M)^Q, 


rf and Ca{r) = T’o ii define a modihed 


h - 2 ^uCoL ■ 


(5.8) 
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Applying the discrete maximal principle to a-Ricci flow (5.4), we can get similar exis¬ 
tence results for constant a-curvature metric. 

Theorem 5.3. Suppose (M, T, 4*) is a weighted triangulated surface. If the initial metric 


r(0) satisfying aRa,i{0) < 0 for all i ^ V, then the normalized a-Ricci flow (5.4) converges 
to a constant a-curvature metric. 


Proof. Note that the a-curvature Ri^a evolves according to (5.5) along the normalized 


a-Ricci flow (5.4). The maximum principle, i.e. Theorem 3.5 is valid for this equation. 


By the maximum principle, if a > 0 and Ra,i < 0 for all z G R, we have 

Rr 


(Ra,min(0) - Ra,av) < Ra,i - Ra,av < Ra,av{'^ “ 


^a,av 


)e 


OiRct^ayt 


-RQ,max (0)' 

a,av _ D \ paRa,avt D . _ U <('R iDi— /? \pO‘Ra,avt 


If a < 0 and Ra,i > 0 for all z G R, we have 
R, 


Ra,min (0) 

In summary, if aRa,i < 0 for all z G R, there exists constants ci and C 2 such that 


cie 


.O^Ra,avi r) , _ p y^^y^^Ra.avt 


^a,av 


< C2e“ 


which implies the exponential convergence of the normalized a-Ricci flow (5.4). 


Using the normalized a-Ricci flow (5.4), the combinatorial a-Calabi flow (5.7) and the 


modified a-Calabi flow (5.8), we can give a characterization of the existence of constant 
a-curvature metric. 

Theorem 5.4. Suppose (M, T, 4>) is a weighted triangulated surface with ax(M) < 0. 
Then the existence of constant a-curvature circle packing metric, the convergence of a- 


Ricci flow (5.4), the convergence of a-Calabi flow (5.7) and the convergence of modified 


a-Calabi flow (5.8) are mutually equivalent. Furthermore, if the solutions of the flows 


converge, then they all converge exponentially fast to a constant R^-curvature metric. 

Proof. Assume there exists a constant curvature metric r*, u* is the corresponding 
u-coordinate. Consider the a-Ricci potential 

N 

(5.9) 


pu ^ 

Fa{u)= '^{Ki - Ra^avrf)dUi. 

i=l 


It has similar properties as stated in Lemma 3.17 Restricted to the hypersurface U = 
{u G R^l Yli we have 

lim Fa{u) = -|-oo 

uGU ,u^oo 
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and Fa\i( is also proper. The rest proof is similar to that of Theorem 3.18 and Theorem 


4.2, so we omit it here. 


We can also use the combinatorial a-Ricci flow, combinatorial a-Calabi flow and mod¬ 
ified combinatorial a-Calabi flow to study the prescribing a-curvature problem. Specifi¬ 
cally, we have the following result. 


Theorem 5.5. Suppose (M, T, 4>) is a weighted triangulated surface, a G lH is a given 
real number, G C{V) is a function defined on M with aRa^i < 0 for i = 1, • • • ,N. 
Then R^ is an admissible a-curvature if and only if the combinatorial a-Ricci flow with 
target 


dvi 

dt 




exists for all time and converges, if and only if the combinatorial a-Calabi flow with target 


= Ao,(Ra - Ra)iri 

exists for all time and converges, if and only if the modified combinatorial a-Calabi flow 
with target 

ii = 

exists for all time and converges, where C = ~ 


Remark 8 . For a = 0, the condition ax(M) < 0 is always satisfied, so there is no 
restriction on x(M) and the equivalences in Theorem 5.4 and Theorem 5.5 are always 
valid. For a = 0, the equivalent condition given by combinatorial Ricci flow is obtained in 
and the equivalent condition given by combinatorial Calabi flow is obtained in 


Now we consider the uniqueness and existence of constant a-curvature metrics. As the 
results are parallel to that of R-curvature, we just give the statements of the results and 
omit the proofs. 


Theorem 5.6. Suppose (M, T, 4>) is a weighted triangulated surface with ax(M) < 0, 
then the constant a-curvature metric is unique if it exists. Specificly, if ax(M) = 0, then 
there exists at most one constant a-curvature metric up to scaling. If ax(M) < 0, then 
for any c*, there exists at most one metric with a-curvature Ra,i = c*. 

Remark 9. For a = 0, this is the result obtained by Thurston [JH] and Chow and Luo 

HU. 
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Theorem 5.7. Suppose is a weighted triangulated surface with ax{M) < 0. 

Then there exists a constant a-curvature metric if and only if there exists a circle packing 
metric r* such that, for any nonempty proper subset I of V, 

*Q; 

27rx(M) ^^y^j >- ^ (7r-<i>(e)) + 27rx(F,). (5.10) 

{e,v)&Lk(l) 

Theorem 5.8. Suppose (M, T, ‘h) is a weighted triangulated surface. If ax{M) < 0, then 
there exists a constant a-curvature metric if and only if n H^q (j) when a > 0 and 
'3^ n / (j) when a < 0. 


Corollary 5.9. Suppose is a weighted triangulated surface with ax{M) > 0. 

If there exists a constant a-curvature metric, then there exists a circle packing metric 
r* such that (5.10) is valid, which implies that 3^ H R^q cj) in the case of a > 0 and 
x{M) > 0, and 3^ n R:^q 4> in the case of a < 0 and x{M) < 0. 


Theorem 5.10. Suppose {M,T, d*) is a weighted triangulated surface. If there is a circle 
packing metric r satisfying > 0 for all i £ V, and 

^ (tt - $(e))-h 27rx(F/) < 0 

{e,v)^Lk(I) 

for any nonempty subset / of R, then there exists a non-negative constant a-curvature 
metric. 


Remark 10. Theorem 3.12 and Theorem 3.20 are now special cases of Theorem 5.3 and 
Theorem 15.101 


6 3-dimensional combinatorial Yamabe problem 

6.1 The definition of combinatorial scalar cnrvature 


Suppose M is a 3-dimensional compact manifold with a triangulation T = {V,E,F,T}, 
where the symbols V,E,F,T represent the sets of vertices, edges, faces and tetrahedrons 
respectively. A sphere packing metric is a map r : R —)• (0, -|-oo) such that the length 
between vertices i and j is kj = ri + rj for each edge {i,j} £ E, and the lengths 
hj,hk, hii Ijk, Ijhlki determines a Euclidean tetrahedron for each tetrahedron {i, j, k, 1} £ T. 
We can take sphere packing metrics as points in R> 0 ) ^ times Cartesian product of (0, oo), 
where N = denotes the number of vertices. It is pointed out |24j that a tetrahedron 
{i,j,k,l} £ T generated by four positive radii ri,rj,rk,ri can be realized as a Euclidean 
tetrahedron if and only if 


Qijkl 



> 0 . 


( 6 . 1 ) 
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Thus the space of admissible Euclidean sphere packing metrics is 

Tlr = {r G I Qijki > 0, V{i, j, k,l} gT } . 

Cooper and Rivin [18] called the tetrahedrons generated in this way conformal and proved 
that a tetrahedron is conformal if and only if there exists a unique sphere tangent to all 
of the edges of the tetrahedron. Moreover, the point of tangency with the edge {i,j} is 
of distance to Vi. They further proved that Tlj- is a simply connected open subset of 
R^q, but not convex. 

For a triangulated 3-manifold (M, T) with sphere packing metric r, there is also the 
notion of combinatorial scalar curvature. Cooper and Rivin [T8| defined the combinatorial 
scalar curvature Ki at a vertex i as angle deficit of solid angles 

Ki = 4:7r- ^ aijki, ( 6 . 2 ) 

where aijki is the solid angle at the vertex i of the Euclidean tetrahedron {i,j,k,l} G T 
and the sum is taken over all tetrahedrons with i as one of its vertices. Ki locally measures 
the difference between the volume growth rate of a small ball centered at vertex Vi in M 
and a Euclidean ball of the same radius. Cooper and Rivin’s definition of combinatorial 
scalar curvature is motivated by the fact that, in the smooth case, the scalar curvature 
at a point p locally measures the difference of the volume growth rate of the geodesic ball 
with center p to the Euclidean ball IMIEI. 

Similar to the two dimensional case. Cooper and Rivin’s definition of combinatorial 
scalar curvature Ki is scaling invariant, which is not so satisfactory. The authors |23] 
once defined a new combinatorial scalar curvature as ^ on 3-dimensional triangulated 

' i 

manifold {M,'T) with sphere packing metric r. Motivated by the analysis in Section]^ we 
find that it is more natural to define the combinatorial scalar curvature in the following 
way. 


Definition 6.1. For a triangulated 3-manifold (M, T) with a sphere packing metric r, 
the combinatorial scalar curvature at the vertex i is defined as 


Ri = 



(6.3) 


where Ki is given by ( 6 . 2 ). 


As Ri differs from Ki only by a factor \, Ri still locally measures the difference between 
the volume growth rate of a small ball centered at vertex Vi in M and a Euclidean ball of 
the same radius. Furthermore, according to the analysis in Section]^ rf is the analogue 
of the smooth Riemannian metric. If = erf for some positive constant c, we have 
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Ri = c~^Ri. This is similar to the transformation of scalar curvature in the smooth case 
under scaling. 

Analogous to the two dimensional case, we can define a measure on the vertices. As 
rf is the analogue of the Riemannian metric, we can take the measure as /U* = r?, which 
corresponds to the volume element. Then the total combinatorial scalar curvature is 

. N N 

s= Rd^i = y^Rirl = y^Kin. (6.4) 

i=i i=i 

Note that S is just the functional introduced by Cooper and Rivin in [18]. For the total 
combinatorial scalar curvature S, we have the following important property. 

Lemma 6.2. ([THl ES] [26]) Suppose {M,T) is a triangulated 3-manifold with sphere 
packing metric r, S is the total combinatorial scalar curvature. Then we have 


= K. 


(6.5) 


If we set 


A = HesSrS 


,Kn) 

d{ri,--- ,rN) 


then A is positive semi-definite with rank N 
spanned by the vector r. 


/ ^ 
' dri 


dKi \ 

dvN 


\ dri 


OKn 
9r-jv / 


1 and the kernel of A is the linear space 


It should be emphasized that , as pointed out by Glickenstein [23], the element 
for i ~ j maybe negative, which is different from that of the two dimensional case. 

By the definition of measure /x, the total measure of M is ^i{M) = We will 

denote the total measure of M by R for simplicity in the following, if there is no confusion. 
The average combinatorial scalar curvature is 


_ Im ^diJ, ^ ^ ^ Kjrj 
- j\,d^ -y- ■ 


( 6 . 6 ) 


6.2 Combinatorial Yamabe problem in 3 dimension 

For the curvature Ri, it is natural to consider the corresponding constant curvature prob¬ 
lem. Suppose Ri = A,Vi G V, for some constant A, then we have Ki = Ar?, which implies 
that X = y. We can take y as a functional of the sphere packing metric r. However, this 
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functional is not scaling invariant in r. So we can modified the functional as yf/a- This 
recalls ns of the smooth Yamabe problem. 

The classical Yamabe problem aims at solving the problem of existence of constant 
scalar cnrvatnre metric on a closed manifold. In order to stndy the constant scalar cur¬ 
vature problem on a closed Riemannian manifold (M,g), Yamabe [50] introduced the 
so-called Yamabe functional Q{g) and the Yamabe invariant Ym,[ 3 o]’ ''^hich are defined as 


Qig) = 


Im ^dVg 


^M,[go] = inf <5(5), 

9^190] 

where [g'o] is the conformal class of Riemannian metric go- Trndinger [39] and Aubin |2] 
made lots of contributions to this problem, and Schoen |36| finally gave the solntion of the 
Yamabe Problem. We refer the readers to |34] for this problem. 

For piecewise flat manifolds, one can introduce similar functionals and invariants. 
Champion, Glickenstein and Young [Hj studied Einstein-Hilbert-Regge functionals and 
related invariants on triangulated manifolds with piecewise linear metrics. The authors 
also introduced a type of combinatorial Yamabe fnnctional and studied its proper¬ 


ties. To study the curvature defined by (6.3), we can introduce the following definition of 


combinatorial Yamabe functional and Yamabe invariant on triangnlated 3-manifolds with 
sphere packing metrics. 


Definition 6.3. Snppose (M, T) is a triangulated 3-manifold with a fixed triangulation 
T- The combinatorial Yamabe functional is defined as 


Q(r) = A = 


r G TIt- 


(6.7) 


The combinatorial Yamabe invariant with respect to T is defined as 


Ym,t = inf Q{r). 

r&XHj- 

The admissible sphere packing metric space TUt- for a given triangulated manifold 
{M,T) is an analogne of the conformal class [g'o] of a Riemannian manifold {M,go), as 
every admissible sphere packing metric could be taken to be conformal to the metric with 
all Tj = 1. We call OJIt- the combinatorial conformal class for (M, 7~). It is uniquely 
determined by the triangulation T of M. Ym,t is referred to as the Yamabe constant for 

(M,r). 
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Note that the combinatorial Yamabe invariant Ym,t is an invariant of the conformal 
class dyij- and is well defined, as we have 


IQ(OI 


5 


T^Kin 

pi/s 


(Sr3)5 



Kh/2, 


( 6 . 8 ) 


and Ki satisfies (4 — 2d)7r < Ki < dvr, where d is the maximal degree and d < E^. If the 
equality in ( |6.8| ) is achieved, the corresponding sphere packing metric must be a constant 
curvature metric. 

By direct computation, we have 

YnQ = — Rav^i ), (6.9) 

which implies that r is a constant combinatorial scalar curvature metric if and only if it 
is a critical point of combinatorial Yamabe functional Q{r). 

Analogue to the smooth Yamabe problem, we can raise the following combinatorial 
Yamabe problem on 3-dimensional triangulated manifold. 


The Combinatorial Yamabe Problem. Given a 3-dimensional manifold M with a 
triangulation T, find a sphere packing metric with constant combinatorial scalar curva¬ 
ture in the combinatorial conformal class 9 JT 7 -. 


We could further consider finding a suitable triangulation for M which admits a con¬ 
stant combinatorial scalar curvature metric. 

It is easy to see that the Platonic solids with tetrahedral cells all admit constant 
combinatorial scalar curvature metric, including the 5-cell, the 16-cell, the 600-cell, etc. 
In these cases, the constant combinatorial scalar curvature metrics arise from symmetry 
and taking the radii equal. 


6.3 Combinatorial Yamabe flow 


In 1980s, Hamilton |30 ( I31| proposed the Yamabe flow to the Yamabe problem. For a closed 
n-dimensional Riemannian manifold {M'^, g) with n >3, the Yamabe flow is defined to be 


with normalization 


m 


II 

1 

(6.10) 

II 

1 

(6.11) 
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where R is the scalar curvature of g and r is the average of the scalar curvature. Along the 
normalized Yamabe flow, the volume is invariant, the total scalar curvature is decreased, 
and the scalar curvature evolves according to 


f) /? 

^ = {n-l)/^R + R{R-r). 


( 6 . 12 ) 


It is proved by Hamilton m that the solution to the Yamabe flow ( |6.11[ ) exists for all time 
and the solution converges exponentially fast to a metric with constant scalar curvature 
if i? < 0 initially. Ye m then proved that, if the initial metric is locally conformally flat. 


then the solution to the normalized Yamabe flow (6.11) converges to a metric of constant 


scalar curvature. Brendle [5] proved that the solution of (6.11) also converges to a metric 
of constant curvature if 3 < n < 5. He [6] further handled the case of n > 6 and got some 
convergence results. 

In the combinatorial case, Lno |35] first introduced the combinatorial Yamabe flow 
on surfaces and Glickenstein |24[ 125] introduced the combinatorial Yamabe flow on 3- 
dimensional manifolds and studied its related properties. Luo [36] also introduced a com¬ 
binatorial curvature flow for piecewise constant curvature metrics on compact triangulated 
3-manifolds with boundary consisting of snrfaces of negative Euler characteristic. To study 
the constant curvature problem of Ri, we introduce the following combinatorial Yamabe 
flow. 

Definition 6.4. Given a 3-dimensional triangulated manifold (M, T) with sphere packing 
metric r, the combinatorial Yamabe flow is defined to be 


with normalization 


dgi _ „ 

dt ~ 


^ = {Ra, - R^)gi, 


(6.13) 

(6.14) 


where gi = rf and Rav is the average of the combinatorial scalar curvature given by (6.6) 


Following the 2-dimensional case, it is easy to check that the solutions of (6.13) and 


(6.14) can be transformed to each other by a scaling procedure. And it is easy to see that. 


if the solution of (6.14) converges to a sphere packing metric r(-|-oo), then r(-|-oo) is a 
metric with constant combinatorial scalar curvature. Analogous to the smooth Yamabe 
flow, we have the following properties. 


Proposition 6.5. Suppose r{t) is a solution of ( |6.14 ) on a triangnlated 3-manifold (M, T). 
Along the flow (6.14), the total measure fi{M) = V = i® invariant and the total 

combinatorial scalar cnrvature S is decreased. 
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Proof. The normalized combinatorial Yamabe flow could be written as 


Using this equation, we have 

dV _ 
dt 


dt ~ 2^^^^ 


= 1^2 - Ri) = 0 . 

i=\ i=l 


For the total combinatorial scalar curvature 5, we have 


dS dKi dri 


dt ^ dt 

2 = 1 


2=1 


N ^ N 

^ ^ ^2 j 2 ^ ^ -^ 2 ) 


1 A Ki{Ki - Rgyrf) 

1 {Ki — Ravf'i)‘^ 

” 2 ^ rv ’ 

2=1 ^ 


2=1 


(6.15) 


which implies that 5 is decreased along the flow (6.14). Note that Lemma 6.2 is used in 
the third step. ■ 


Remark 11. As the total measure of M is invariant along (6.14), we will focus on the 
properties of (6.14) in the following. Furthermore, we will assume r(0) G = {r G 

Iklls = ™ following. 

By direct calculations, we find that the curvature Ri evolves according to the following 
equation along the normalized combinatorial Yamabe flow (6.14) 


dR, 


1 Ask 




2rf ^ drj 

2 j = l J 


(6.16) 


If we define the Laplacian as 


1 Ask 


= -3 EA = 3 E(-|33)(/, - /.) 


r; ^ Sr,- 

I j 


J~2 


(6.17) 
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for / G C{V), then the equation (6.16) could be written as 


^ = ]^i^Ri + Ri{Ri-Ray), 


(6.18) 


which has almost the same form as the evolution equation (6.12) of scalar curvature along 
the Yamabe flow (6.11) on three dimensional smooth manifolds. The Laplacian defined 
by (6.17) satishes A/d/r = 0 and Ac = 0 for any / G C{V) and constant c G H. 
Furthermore, by the calculations in [2l], we have 


• j ■ ■ 


(6.19) 


where l*j is the area dual to the edge {i,j}- Note that this form of Laplacian is very 


similar to Hirani’s definition |32j of Laplace-Beltrami operator 


A/i = 


1 


I*. 


V* ^ I 






except the first factor. It should be mentioned that rf is a type of volume. 

Though the definition (6.17) of Laplacian has lots of good properties, it is not a Lapla¬ 
cian on graphs in the usual sense, as l*j may be negative, which makes the maximum 
principle for (6.18) not so good. This was founded by Glickenstein and he studied the 
properties of such Laplacian in [JU |25] . 


To study the long time behavior of (6.14), we need to classify the solutions of the flow. 


Definition 6.6. A solution r(t) of the combinatorial Yamabe flow (6.14) is nonsingular 
if the solution r(t) exists for t G [0, -|-oo) and {r(t)} CC Tlj- n 

In fact, the condition {r(t)} CC 9 II 7 - Pi ensures the long time existence of the 

flow (6.14). Furthermore, we have the following property for nonsingular solutions. 


Theorem 6.7. If there exists a nonsingular solution for the flow (6.14), then there exists 
at least one sphere packing metric with constant combinatorial scalar curvature Ri on 


Proof. By (6.15), we know that Q{r) is decreasing along the flow (6.14). As Q{r) 
is uniformly bounded by (6.8), the limit hmt_,._|_oo (5(r’(t)) exists. Then there exists a 
sequence tn t +00 such that 


N 




„ 2\2 


i=l 


Vi 


( 6 . 20 ) 
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As {r’(i)} CC OJlr n there exists a subsequence, denoted as r„, of r{tn) such that 

rn ^ r* e OJIt- bl Then (6.20) implies that r* satisfies K* = Rav{r*)'^ and r* is a 

sphere packing metric with constant combinatorial scalar curvature Ri. ■ 


In fact, under some conditions, we find that the sphere packing metrics with constant 
combinatorial scalar curvature are isolated on 

Theorem 6.8. The sphere packing metrics with nonpositive constant combinatorial scalar 
curvature are isolated in SOIt- n 

Proof. Define 

r ^ (—{Ki - Xrl), • • • , —{Kn - Xr%) \ , 

Vn tn ) 

where A = ^. It is easy to see that the zero point of G corresponds to the metric with 
constant combinatorial scalar curvature. By direct calculations, the Jacobian of G at the 
constant curvature metric is 

DG= - 2ASi E-i, (6.21) 

where If r* is a sphere packing metric with nonpositive constant 

combinatorial scalar curvature, then A < 0 and hence DG is a positive semi-definite matrix 
with rank N — 1 and kernel {cr^|c G R}, which is normal to Restricted to 

DG is positive definite and then nondegenerate, which implies that the zero points of G 
with nonpositive curvature are isolated in dJl-j- n ■ 


Suppose the solution of the flow (6.14) exists for t G [0, T) with 0 < T < + 00 , then we 
have {r(t)} C SOIt- n However, if 


{r(t)}n9(9JTrn§^-^) ^0, 


the flow will raise singularities. They could be separated into two types as follows. 

Essential Singularity There is a vertex i G R such that there exists a sequence of time 
tn\ T such that ri{tn) —t 0 as n —>■ + 00 ; 

Removable Singularity There exists a sequence of time tn ^ T such that ri{tn) CC 
R>o for all vertices i G R, but there exists a tetrahedron {ijkl} G T such that 
Qijkiitn) -t 0 as n + 00 . 
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It seems that the following conjectures are likely to hold for the normalized combina¬ 
torial Yamabe flow (6.14) on 3-dimensional manifolds. 


Conjecture 1. The normalized combinatorial Yamabe flow (6.14) will not develop essen¬ 
tial singularity in finite time. 


Conjecture 2. If no singularity develops along the normalized flow (6.14), the solution 


converges to a sphere packing metric with constant combinatorial scalar curvature as time 
approaches inhnity. 


It is interesting to note that Glickenstein [25] made a small amount of progress on 
Conjecture 1. He proved that, for some special class of complexes, the flow he introduced 
would not develop singularities in finite time. 

We have mentioned above that the existence of sphere packing metric with constant 
combinatorial scalar curvature is necessary for the convergence of the normalized com¬ 


binatorial Yamabe flow (6.14). In fact, it is almost sufficient and we have the following 
result. 

Theorem 6.9. Suppose r* is a sphere packing metric on (M, T) with nonpositive con¬ 
stant combinatorial scalar curvature. If ||r(0) — r*|p is small enough, the solution of the 
normalized combinatorial Yamabe flow (6.14|) exists for all time and converges to r*. 


Proof. We can rewrite the normalized combinatorial Yamabe flow (6.14) as 


dr. 


1 


dt - 

Set rj(r) = \{Rav — Ri)ri, then by the calculations in the proof of Theorem 


6.8 


we have 


/ — 


DT\r* = — 


rjv 


where 


Z = A-2Ra 


r'J^ 




1/2 


I - 


^ 3/2 . (-^ 3 / 2^1 

V 






1/2 
N ) 


If Rav < 0, then DT\r* is a matrix with rank N — 1 and kernel {cr^|c G IR}. Furthermore, 
the nonzero eigenvalues of DT\j.* are all negative. Note that, along the normalized flow 


(6.14), ||r ||3 is invariant and thus the kernel {cr^|c G K,} is transversal to the flow. This 
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implies that DT\r* is negative definite on ^ and r* is a local attractor of the normal¬ 


ized combinatorial Yamabe flow (6.14). Then the conclusion follows from the Lyapunov 
Stability Theorem. ■ 


This theorem has the following interesting corollary. 

Corollary 6.10. Given a 3-dimensional triangulated manifold (M^, T). If the initial total 
combinatorial scalar curvature functional 5(0) < 0, and no singularity develops along the 


normalized combinatorial Yamabe flow (6.14), then the solution converges to a sphere 


packing metric r* with nonpositive constant combinatorial scalar curvature. 


Proof. By the proof of Theorem 6.7 if the solution r{t) of (6.14) is nonsingular, there 
exists a sequence t r*, where r* is a sphere packing metric with constant combinatorial 
scalar curvature. As S is decreasing along the flow ( |6.14[ ) and 5(0) < 0, r* has nonpositive 
combinatorial scalar curvature. Then Theorem |6.9| implies the conclusion of the corollary. 


Remark 12. Especially, if 72* (0) < 0, Corollary 6.10 is still valid. 


7 Some questions 

There are several questions which we find interesting relating to the results in the paper. 

1. From the results in Section we see that the surfaces with positive Euler char¬ 
acteristic are particular for the constant curvature problem on surfaces. Example 
shows that combinatorial Ricci flow (1.12) may not converge to a constant 72- 
curvature metric, while Example shows that constant curvature metric may not 
be unique. If we want to approximate a smooth geometric object by corresponding 
discrete objects, the quantities on these objects should obey similar laws or exhibit 
similar properties. However, it has been proved [sniES] that the Ricci flow on sur¬ 
faces converges to a constant curvature metric for any initial metric. The differences 
between discrete case and the smooth case are very interesting. It deserves deeper 
studying. Maybe these differences are caused by the triangulation. The tetrahedron 
triangulation used in Example and maybe too rough to approximate the sphere. 
We expect these differences will disappear as the triangulation becomes finer and 
finer. We want to know whether we can find a triangulation for the sphere such 
that the constant curvature metric uniquely exists. We further want to know how 
to evolve discrete curvature 72* = 7L*/r? along discrete curvature flows to constant 
curvature or, more generally, evolve discrete a-curvature to constant curvature when 
axiM) > 0 . 
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2. Whether there are similar topological and combinatorial obstructions, similar to 
Thurston’s criterion (3.22), for the existence of constant a-curvature metric? For 


the negative constant curvature metric, we derived an existence result (Theorem 


3.12) by the discrete maximum principe, which is different from the other results. 


We want to further develop these methods to derive more existence results. We want 
to know whether we can get similar convergence results as that of the smooth surface 
Ricci flow and then give a discrete uniformization theorem for Rj. We believe that 
the results in [1] will play an important role in the procedure if it is feasible. We 
also want to know the relationship between our conditions and Thurston’s criterion 


(3.22). 


3. The first positive eigenvalue of discrete Laplace operator plays an important role in 
the proof of the main results in two dimension. It is interesting to estimate the first 
positive eigenvalue of discrete Laplace operator and use it to derive some existence 
results. 

4. Investigate the prescribing curvature problem when the target curvature Ri > 0 or, 
more generally, aRa,i > 0 for at least one vertex i. 


5. Study the singularities of three dimensional combinatorial Yamabe flow (I.I9), and 


find topological and combinatorial obstructions for the existence of sphere packing 
metric with constant combinatorial scalar curvature. 


6. Study the discrete Yamabe functional (6.3) and solve the combinatorial Yamabe 
problem in three dimension. 

7. Study the combinatorial Yamebe problem in higher dimensions. 
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